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[i] This study presents a global assessment of the sensitivity of droplet number to 
diabatic activation (i.e., including effects from entrainment of dry air) and its first-order 
tendency on indirect forcing and autoconversion. Simulations were carried out with the 
NASA Global Modeling Initiative (GMI) atmospheric and transport model using 
climatological metereorological fields derived from the former NASA Data 
Assimilation Office (DAO), the NASA Finite volume GCM (FVGCM) and the 
Goddard Institute for Space Studies version II’ (GISS) GCM. Cloud droplet number 
concentration (CDNC) is calculated using a physically based prognostic 
parameterization that explicitly includes entrainment effects on droplet formation. 

Diabatic activation results in lower CDNC, compared to adiabatic treatment of the 
process. The largest decrease in CDNC (by up to 75%) was found in the tropics and in 
zones of moderate CCN concentration. This leads to a global mean effective radius 
increase between 0.2-0.5 jrm (up to 3.5 //m over the tropics), a global mean 
autoconversion rate increase by a factor of 1.1 to 1.7 (up to a factor of 4 in the tropics), 
and a 0.2-0.4 W m 2 decrease in indirect forcing. The spatial patterns of entrainment 
effects on droplet activation tend to reduce biases in effective radius (particularly in the 
tropics) when compared to satellite retrievals. Considering the diabatic nature of ambient 
clouds, entrainment effects on CDNC need to be considered in GCM studies of the 
aerosol indirect effect. 

Citation: Barahona. D., R. Sotiropoulou, and A. Nenes (2011), Global distribution of cloud droplet number concentration, 
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1. Introduction 

[ 2 ] Cloud droplet number concentration (CDNC) depends 
on the size distribution and composition of the precursor 
aerosol, and, on the thermodynamical and dynamical (i.e., 
updraft velocity, mixing rates) state of the cloudy air during 
its formation [Seinfeld and Pandis, 1998]. Increasing the 
precursor aerosol number concentration may lead to a 
decrease in cloud effective radius and therefore to an 
increase in cloud albedo, i.e., the first aerosol indirect effect 
[Twomey, 1977]. The limited ability of general circulation 
models (GCMs) to explicitly represent aerosol-cloud inter- 
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actions motivated the development of physically based 
representations parameterizations of the cloud droplet 
formation process [Abdul -Razzak and Ghan, 2000; Nenes 
and Seinfeld, 2003; Fountoukis and Nenes, 2005; Ming et al., 
2006; Kumar et al., 2009]. State-of-the-art parameteriza- 
tions take into account effects of the aerosol composition, 
such as the role of organic surfactants and solutes [Fountoukis 
and Nenes, 2005; Abdul -Razzak and Ghan, 2005], adsorption 
activation [Kumar et al, 2009], and kinetic and mass transfer 
limitations on droplet formation [Fountoukis and Nenes, 2005; 
Asa-Awuku and Nenes, 2007; Barahona et al., 2010], and 
mixing and entrainment [Barahona and Nenes, 2007]. 

[ 3 ] State-of-the-art atmospheric models use double- 
moment cloud microphysics schemes (which vary signifi- 
cantly in how each cloud moment is computed), in which the 
time evolution of the moments of the droplet size distribution 
is determined by grid-scale advection, detrainment, turbulent 
diffusion, and microphysical processes [e.g., Takemura 
et al., 2005; Quaas et al., 2009; Chen et al., 2010; Ghan 
et al., 1997; Lohmann and Feichter, 1997; Morrison and 
Gettelman, 2008]. The calculated total liquid condensate 
mass (LWC) and the cloud droplet number concentration 
(CDNC) are then used to estimate the cloud optical properties 
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(e.g., droplet effective radius, R e fj, and optical depth) and 
autoconversion rate. In these schemes cloud condensation 
nuclei (CCN) activation represents a source of cloud droplets. 
CCN activation is typically parameterized considering the 
adiabatic ascent of air parcels from cloud base [e.g., Abdul - 
Razzak and Ghan, 2000; Fountoukis and Nenes, 2005; 
Ming et al., 2006; Segal and Khain, 2006]. It is, however, 
well known that real clouds are predominantly subadiabatic; 
field observations show both CDNC and LWC below the 
adiabatic limit, which results from mixing and entrainment of 
dry air throughout the cloud lifetime [Brenguier and Chaumat, 
2001; Conant et al., 2004; Lu et al., 2008; Pruppacher and 
Klett, 1997], 

[ 4 ] Mixing and entrainment of environmental air affects 
cloud formation over multiple spatial scales. Radiatively 
induced turbulent mixing near cloud top may trigger the 
breakup and dissipation of stratocumulus decks [e.g., 
Bretherton et al., 2004; Stevens, 2005; Gerber et al., 2005], 
Lateral entrainment weakens updrafts affecting the forma- 
tion of cumulus and deep convective clouds [e.g., Carpenter 
et al., 1998; Chosson et al., 2007; Cohen, 2000; Lin and 
Arakawa, 1997; Rogers et al., 1985; Stommel, 1947]. The 
microphysics of cloud formation is also affected by 
entrainment as it modifies the drop size distribution [e.g., 
Pruppacher and Lee, 1976; Paluch, \919\Ragaetal, 1990; 
Raga and Jonas, 1993; Cohen, 2000; Neggers et al., 2003; 
Barahona and Nenes, 2007; Lu et al., 2008]; neglecting 
entrainment effects on CDNC (the current assumption in 
GCMs) produces much more narrow droplet size dis- 
tributions when compared to observations [e.g. Lin and 
Arakawa, 1997; Cohen, 2000; Brenguier and Chaumat, 
2001; Snider et al., 2003; Morrison and Grabowski, 
2008; Pruppacher and Klett, 1997; Hsieh et al, 2009b]. 

[ 5 ] GCM cloud schemes account implicitly for entrain- 
ment and mixing effects on liquid condensate by introduc- 
ing empirically derived subgrid probability distributions of 
LWC or by prescribing the LWC vertical profile as a 
function of the large-scale temperature and relative humidity 
fields [e.g., Takemura et al., 2005; Quaas et al., 2009; Chen 
et al., 2010; Ghan et al., 1997; Lohmann and Feichter, 
1997; Morrison and Gettelman, 2008; Flack, 1998]. How- 
ever, the CDNC vertical profile is often still assumed equal 
to the upper limit of adiabatic activation (Figure lb). This 
treatment is inherently inconsistent, with a tendency for over- 
estimating CDNC and underestimating droplet size, with 
important implications for indirect effect assessments [e.g., 
Kim et al. 2008], 

[6] The assumption of adiabatic activation stems from 
considering droplet formation in large-scale clouds as 
occurring solely within the core of horizontally homoge- 
neous stratocumulus decks; this approach, however, leaves 
out the very common trade wind cumulus of the tropical 
regions [Albrecht et al., 1995]. Furthermore, it is known that 
in tropical and subtropical marine environments cumulus 
originating at the surface mixed layer are responsible for the 
maintenance of stratocumulus by supplying liquid water to 
the cloud layer [ Albrecht et al., 1995; Bretherton and 
Pincus, 1995], Stratocumulus are known to cycle between 
horizontally inhomogeneous and homogeneous states, the 
former resembling trade wind cumulus fields [Albrecht et al., 
1995]. These transitions can be induced by cumulus forming 
below the cloud layer, cloud top entrainment, and drizzle 


[Albrecht et al., 1995; Feingold et al., 2010; Korolev, 1995; 
Bretherton and Pincus, 1995], Cloud droplet formation in 
shallow cumulus is more strongly influenced by entrainment 
than in the core of stratocumulus decks [Pruppacher and 
Klett, 1997; R. Morales et al., Evaluation of a diabatic 
droplet activation parameterization using in-situ cloud data, 
submitted to Journal of Geophysical Research, 2010]. Thus, 
in real clouds, diabatic cloud droplet formation occurs much 
more frequently than adiabatic [e.g., Chosson et al., 2007; 
Gerber et al., 2005; Grabowski and Pawlowska, 1993; Lin 
and Arakawa, 1997; Wang and Albrecht, 1994; Morales 
et al., submitted manuscript, 2010], which is not recognized 
in most GCM studies of the indirect effect. 

[ 7 ] Although there is a considerable body of work on 
entrainment effects on large-scale cumulus convection 
parameterization [e.g., Morton et al, 1956; Kain and Fritsch, 
1990; Siebesma and Cuijpers, 1995; Randall et al., 2003; 
Donner et al., 2001] less has been done on linking mixing 
processes with the microphysical parameters (primarily 
CDNC) of clouds in GCMs. Ideally, one would like to 
describe all of possible entrainment processes during the 
cloud lifetime (e.g., homogeneous [Pruppacher and Klett, 
1997] and inhomogeneous mixing across all scales, and the 
droplet distribution changes that result [e.g., Pruppacher and 
Klett, 1997; Brenguier and Chaumat, 2001; Morrison and 
Grabowski, 2008]) (Figure la). This is very difficult, how- 
ever, to determine, especially in a prescribed manner for 
inclusion in a subgrid parameterization. There may be a way, 
however, to describe the average cloud column properties that 
arise from the effects of cloud base activation and subsequent 
dilution and drying of droplets from entrainment. Toward 
this, Barahona and Nenes [2007] noted that mixing during 
the droplet activation process acts to decrease the cloud 
maximum supersaturation (hence activated CDNC) com- 
pared to adiabatic activation (Figure lc). Adopting a homo- 
geneous entrainment and mixing cloud model, Barahona and 
Nenes [2007] developed a droplet formation parameteriza- 
tion, using an “effective” per length entrainment rate, e that 
accounts for the effects of nucleation near cloud base and 
deactivation further up in the cloud column (Figure lc). 
Using scaling arguments and microphysical observations of 
nonprecipitating cumulus, Morales et al. (submitted manu- 
script, 2010) concluded that the [Barahona and Nenes, 
2007] formulation can describe cloud-average CDNC in 
entraining cumulus clouds when the average adiabaticity 
ratio (defined as the ratio of cloud liquid water content over 
the adiabatic value) is used to diagnose e. The same study 
also showed that for thick stratocumulus clouds, either 
adiabatic or diabatic formulations predict pdf-average CDNC 
to within experimental accuracy. Therefore, shallow and 
trade wind cumulus and stratocumulus are the main cloud 
types modified by this approach, the effects of which are 
most prominent in the tropics. 

[8] In this work, we study the potential impact of 
entrainment on global CDNC, autoconversion rate, effective 
radius, and the first aerosol indirect effect [Twomey, 1977]. 
This is accomplished by implementing the Barahona and 
Nenes [2007] parameterization in the NASA Global Mod- 
eling Initiative (GMI) chemistry and transport model 
[Meskhidze et al, 2007]. The sensitivity of global CDNC, 
cloud effective radius (R e ff), indirect forcing (IF) and auto- 
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Figure 1. CDNC and LWC calculation in GCMs. (a) Representation of the entrainment profile in an 
ambient cloud where both LWC and cloud droplet activation are affected by entrainment, (b) Current 
GCM representation assuming adiabatic activation but LWC affected by entrainment, (c) The activation 
approach proposed in this study, where e is such so the cloud-average CDNC is similar to the ambient 
cloud in Figure la. Grey arrows represent entrainment of dry air. 


conversion to diabatic activation is studied using different 
meteorological fields to drive the GMI model. 


2. Model Description 

[ 9 ] The NASA GMI model is a highly modular 3-D 
chemistry and transport model [Rotman et al., 2001] capable 
of multiyear simulations for assessments of anthropogenic 
impacts on atmospheric chemistry. GMI computes the first- 
order response (i.e., without feedbacks of the cloud micro- 
physics on the large-scale fields) of climate to external 
forcing and is therefore suitable for studying the first aerosol 
indirect effect [Twomey, 1977]. The modeling framework 
used is described by Meskhidze et al. [2007] and is summa- 
rized here. The aerosol module is coupled to the GMI-CTM 


advection core. Aerosol microphysics is not considered. 
Externally mixed prognostic aerosol species include sulfate, 
black carbon and organic matter from biomass burning and 
fossil fuel, organic matter from natural sources, mineral dust, 
and sea salt [Liu et al., 2005, 2007]. The model considers 
primary emissions, chemical production of sulfate, gravita- 
tional sedimentation, dry deposition, wet scavenging, and 
hygroscopic growth [Liu et al., 2005], Anthropogenic and 
natural aerosol precursor emissions include SO 2 , organic 
matter, black carbon, oceanic DMS, dust and sea salt [Liu 
et al., 2005; Meskhidze et al., 2007]. 

[ 10 ] Two types of clouds exist in GMI: convective and 
large-scale (stratocumulus and trade wind cumulus); clouds 
are allowed to form in any model layer with the exception of 
the layer nearest to the surface, as it is very shallow (formation 
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of fog is not considered). Large-scale cloud fraction is 
diagnosed based on the grid cell relative humidity [Hack, 
1998], RH, using a threshold relative humidity for con- 
densation specified as a function of pressure [Xu and 
Krueger, 1991], For convective cloud fraction, a parame- 
terization with convective mass flux [Xu and Krueger, 
1991] was adopted. Total cloud fraction in each layer is 
obtained from the combination of the large-scale cloud 
fraction and convective cloud fraction assuming random 
mixing [Feng et al., 2004], which has been shown to 
adequately reproduce the observed global distribution of 
liquid water path. Following Hack [1998], the cloud liquid 
water content (LWC) is vertically distributed using a cloud 
liquid water density profile [Kiehl, 1994] derived from 
experimental measurements and GCM simulations. In-cloud 
liquid water path (LWP) at each level is found by the inte- 
gration of LWC between layers [Meskhidze et al., 2007]. 
Cloud forcing and radiative properties are calculated online 
using the CLIRAD-SW solar radiative transfer model [Chou 
et al., 1998], which gives the cloud optical depth and the 
shortwave radiative flux from the surface to the top of the 
atmosphere (TOA). Radiative fluxes are integrated over 
the solar spectrum, from 0.175 pm to 10 pm. Inputs in 
the radiative transfer module include the cloud droplet 
effective radius, the in-cloud LWC, the cloud fraction, the 
specific humidity and the 0 3 concentrations for each grid 
box as well as the surface albedo for the direct and dif- 
fuse light and the solar zenith angle. Aerosol indirect 
forcing (IF) is calculated as the difference in the TOA net 
outgoing shortwave radiative flux between two simula- 
tions, usually between those with “present-day” (natural 
and anthropogenic) and “preindustrial” (only natural) 
emissions of aerosols (and their precursors). In computing 
indirect forcing, surface albedo for direct and diffuse light 
was obtained from the GEOS4 GCM [Bloom et al., 
2005]. Following Del Genio et al. [1996], liquid clouds 
were those with a cloud top temperature above 263.15 K 
over land, and, 269.15 K over ocean. The sensitivity of 
indirect forcing to the meteorological field (used to 
diagnose grid cell RH, T and p) is assessed by carrying 
out simulations with fields from the former NASA Data 
Assimilation Office (DAO), the NASA finite volume 
GCM (FVGCM), and the Goddard Institute for Space 
Studies version II’ (GISS) GCM. Each of the archived 
meteorological data sets spans over 1 year. 

[n] The horizontal resolution in all simulations is, 
4° latitude by 5° longitude. Vertical resolution differs in 
each meteorological field, being 46, 42, and 23 layers for 
DAO, FVGCM, and GISS, respectively. Monthly and 
annually averaged CDNC, R e g, and autoconversion rates 
(Q alU ) are computed from the simulations. Simulations with 
GMI are carried out for each meteorological field (DAO, 
GISS, FVGCM) using both “present day” (i.e., all emissions 
active) and “preindustrial” (i.e., emissions active except 
anthropogenic) emissions. 

2.1. Description of Cloud Parameterizations 

[ 12 ] CDNC is assumed equal to the nucleated droplet 
number which is supported by field campaign studies [e.g., 
Conant et al., 2004; Guibert et al., 2003; Lu et al., 2008; 
Meskhidze et al., 2005; Peng et al., 2005] and justified as 
precipitating cumulus and stratocumulus clouds typically 


have short lifetimes [Pruppacher and Klett, 1997]. Cloud 
droplet number is calculated using the parameterization 
of Barahona and Nenes [2007, hereafter BN07], which 
explicitly includes the effects of entrainment and mixing on 
maximum supersaturation and CDNC. BN07 is an extension 
of the works of Nenes and Seinfeld [2003] and Fountoukis 
and Nenes [2005]; it can treat the effects of externally 
mixed aerosol, CCN containing surfactants and slowly 
growing droplets (expressed by changes in the water vapor 
mass transfer coefficient). This parameterization has been 
shown to reproduce in situ measurements of CDNC in 
cumulus and stratocumulus clouds [Meskhidze et al, 2005; 
Fountoukis et al., 2007; Morales et al., submitted manu- 
script, 2010]. In this study, the effect of aerosol equilibrium 
liquid water content on supersaturation is neglected (i.e., the 
“simplified” version of BN07 is used). This assumption 
introduces less than 1% error in CDNC calculation for 
pristine and moderately polluted environments (aerosol 
mass less than 102 /ig m~ 3 ) and cloud-free air conditions far 
from saturation (RH < 95%), which characterizes most of 
the atmosphere. With these assumptions, BN07 involves 
finding the root of the local supersaturation tendency 
equation [Barahona and Nenes, 2007], 




(1) 


where V is the updraft velocity, s is the supersaturation, ^ is 
the rate of condensation of liquid water onto the drops, 
determined using the method of “population splitting” as 
described by Nenes and Seinfeld [2003], Fountoukis and 
Nenes [2005], and Barahona and Nenes [2007], The criti- 
cal entrainment rate, e c , is defined as the entrainment rate 
that completely prevents droplet activation (i.e., the upper 
limit in e for which cloud formation is possible) given by 
[Barahona and Nenes, 2007], 




a 


(1 -RH) 


AH„M W AT 
RT - 


( 2 ) 


where a = ^ , RH is the ambient relative humidity 

and AT is the difference between cloud and ambient tem- 
peratures. Parcel model simulations show that entrainment 
would affect CCN activation if e/e c >0.1 [Barahona and 
Nenes, 2007]; adiabatic activation occurs when e = 0 (i.e., 
negligible mixing rate). The value of s determined from 
equation (1) corresponds to the maximum supersaturation, 
s max . The droplet number concentration is then calculated as 
the CCN with critical supersaturation less than s max . BN07 
evaluated the parameterization against a detailed numerical 
parcel model over a wide set of cloud forming conditions 
that included entrainment; the mean relative error in com- 
puted CDNC was 2 ± 21%. 

[ 13 ] In principle, equation (1) is also applicable to 
homogeneous mixing in convective clouds, and may pro- 
vide a simple way to parameterize entrainment in large-scale 
convection schemes. Equation (1) is derived assuming that 
mixing between cloudy and cloud-free air is instantaneous 
and homogeneous. In our approach, e and e/e c therefore 
represent “effective” mixing rates (Figure 1) that would 
account for the net effect of mixing processes throughout the 
cloud column, i.e., the reduction in CDNC with respect to 
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Table 1. Dry Aerosol Size Distributions Used in This Study [Lance et al., 2004] a 


Aerosol 

Mode 


Marine 



Continental 


Nuclei 

Accumulation 

Coarse 

Nuclei 

Accumulation 

Coarse 

N t 

0.561 

0.432 

0.008 

0.555 

0.444 

4 x 10~ 4 

4. 

0.02 

0.092 

0.58 

0.01 

0.067 

0.93 

<7 g i 

1.47 

1.6 

2.49 

1.6 

2.1 

2.2 

^ sol 

0.33 

0.33 

0.95 

0.5 

0.5 

0.5 


a D gi , a gh and c so i is the geometric mean diameter (/an), spectral width, and soluble fraction, respectively. A,- is the number concentration of mode “i”, and 
N lo , is total number concentration in the population. 


adiabatic activation. This is justified as the GMI model si- 
mulates cloud-scale averages over relatively large grid cells 
and time steps. From the CDNC perspective, this approach 
captures cloud-column average CDNC better than the adi- 
abatic value in trade wind cumulus and stratocumulus 
(Morales et al., submitted manuscript, 2010). The extension 
of this approach when inhomogeneous mixing and partial 
droplet evaporation significantly impact CDNC requires an 
analysis using ambient cloud data and will be the subject of 
a future study. 

[ 14 ] Autoconversion rate of cloud water to rain is com- 
puted from the scheme of Khairoutdinov and Kogan [2000], 
which is derived from large eddy simulations of drizzling 
stratocumulus clouds. Autoconversion rate (s^ 1 ) is given by 

Qaut = 1350 (J^ (CDNC)- 1 - 79 (3) 

where Q is the cloud fraction and qi the liquid water mass 
mixing ratio; CDNC is in cm~ 3 . Hsieh et al. [2009a] showed 
that among several parameterizations commonly used in 
GCM studies, equation (3) has the strongest dependency on 
CDNC, hence it can express the upper limit of auto- 
conversion sensitivity to CDNC from entrainment. Q aut is 
used as a diagnostic of the potential effect of diabatic acti- 
vation on the cloud lifetime and precipitation rate. However, 
as no feedback is considered between CDNC and LWC, 
Q aut is only meaningful for the comparison between diabatic 
and adiabatic cases. 

2.2. Implementation of the Cloud Activation 
Parameterization 

[ 15 ] Calculation of CDNC using BN07 requires the 
knowledge of T, p, aerosol size distribution and composi- 
tion, V, e, AT and cloud-free air RH. The characteristic 
cloud formation T and p [Barahona and Nenes, 2007] are 
taken as the grid-cell average T and p. Grid-cell average RH 
is assumed to be representative of the cloud-free air en- 
trained in the cloudy columns. AT is estimated from the 
release of latent heat during condensation [Paluch, 1979; 
Seinfeld and Pandis, 1998]; an energy balance applied to a 
cloudy parcel at steady state then gives, 



Application of equation (4) provides a global average AT ~ 
0.4 K for the three meteorological fields (DAO, G1SS, and 
FVGCM), with AT in marine environments being slightly 
higher (~0.5 K) due to the larger LWC in these regions (not 
shown). 


[ 16 ] Calculation of CDNC requires knowledge of the 
aerosol size distribution and hygroscopicity. As currently 
GMI simulates only aerosol mass, distributions are diag- 
nosed from the online simulation by scaling distributions 
obtained from observations [Lance et al., 2004] using the 
predicted sulfate mass. In these distributions, aerosol over 
marine regions is composed of 67% insoluble material and 
33% ammonium sulfate in the fine mode, and, 95% sea salt 
and 5% insoluble material in the coarse mode. Aerosol over 
land is composed of 50% insoluble material and 50% 
ammonium sulfate in both the coarse and fine modes [Lance 
et al, 2004]. The size distribution mode diameter, geometric 
standard deviation, and number concentration for marine and 
continental aerosol are given in Table 1. Chemical effects on 
hygroscopicity are also important [Lance et al., 2004; Nenes 
et al., 2002; McFiggans et al., 2005]; however, they are not 
considered in this study to facilitate the interpretation of the 
effects of entrainment on droplet activation. The effect of 
assuming a constant aerosol size distribution is discussed in 
section 4. 

[ 17 ] This study adopts two approaches to diagnose e. In 
the first one, adiabaticity (i.e., the ratio of in-cloud LWC to 
its adiabatic value; [Kim et al., 2008]) is assumed constant 
and used as proxy for e/e c . Morales et al. (submitted man- 
uscript, 2010) have shown that this approach can give good 
CDNC closure if 1 - e/e c is set equal to the cloud-average 
adiabaticity. Recent in situ observations report -0.5-0.6 
adiabaticity for small cumulus clouds [Lu et al., 2008; 
Morales et al., submitted manuscript, 2010]; thus entrain- 
ment effects on CDNC, IF, and Q aut are studied by carrying 
out simulations with effective e equal to 0.0e c (adiabatic), 
0.4e c , and 0.6e c . We also directly specify e, using average 
values derived from LES simulations of shallow convective 
clouds [Siebesma and Cuijpers, 1995]; e = 2 x 1(T 3 m _1 for 
positively buoyant cloudy parcels (“core” case), and, e = 3 x 
10 3 m 1 for all cloud parcels with positive LWC and 
V (“updraft” case). Similar e values have been reported in 
other published studies that reflect changes in a variety of 
dynamical conditions [e.g., McCarthy, 1974; Raga et al., 
1990; Gregory, 2001; Neggers et al., 2003], so e — 2 x 
1CT 3 - 3 x 1CT 3 nT 1 can be considered representative of the 
global average. Using a fixed, average e/e c implies a spatial 
variation in e (as e c depends on local RH and T, equation (2)) 
whereas maintaining e constant implies a spatial variation in 
e/e c \ thus, the two approaches used represent limits of vari- 
ability in both e and e/e c . 

[is] In addition to the aerosol size distribution and e, 
CDNC calculation requires the knowledge of cloud-scale V. 
Since this is not directly available in large-scale models, we 
prescribe V using values from field campaign studies; an 
average V is thus prescribed for “continental” and “marine” 
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Table 2. Global Mean CDNC, R e jj, Q aut , and IF for Runs Using 
Adiabatic Activation a 



DAO 

GISS 

FVGCM 

< 7 

CDNC (cm 3 ) 

102.6 

96.3 

95.6 

21.4 

Reff (pm) 

8.2 

8.9 

7.5 

0.37 

log (8 au t) 

-10.3 

-10.2 

-10.4 

0.32 

IF (W m 2 ) 

-1.28 

-1.30 

-1.75 

0.33 


a The standard error, a, was calculated as the global mean standard 
deviation from the average of the three meteorological fields. 


clouds equal to 1 m s _1 and 0.35 m s -1 , respectively [Lance 
et al., 2004; Peng et al., 2005; Fountoukis et al, 2007] for 
droplet activation in stratocumulus. A single V can express 
the vertical velocity distribution average droplet number 
in cumulus and stratocumulus clouds, provided that it 
expresses the average updraft velocity distribution in the 
boundary layer [Lance et al., 2004; Meskhidze et al., 2005; 
Peng et al., 2005; Fountoukis et al., 2007; Morales and 
Nenes, 2010]. Although this approach can provide realistic 
droplet distributions [Chen et al., 2010]; vertical velocity is 
nevertheless treated as an adjustable parameter. 


3. Results 

[ 19 ] Entrainment reduces s max hence the number of acti- 
vated CCN (CDNC). The extent of the CDNC reduction 
(compared to assuming adiabatic activation) depends pri- 
marily on e/e c and the characteristics of the aerosol popu- 
lation [Barahona and Nenes, 2007]. Simulations are carried 
out varying e (using both methods outlines in section 2.2) in 
the CDNC parameterization but maintaining the (diabatic) 
LWC as given by GM1 (section 2). This allows quantifying 
the error associated with neglecting entrainment effects on 
CDNC, R ef/ , IF and Q au , 

[ 20 ] Entrainment effects are expressed in terms of an 
absolute change, 

AX ahs =X-X ad (5) 

where X refers to the annual average of the variable being 
analyzed (e.g., CDNC, IF) calculated for diabatic activation, 
and X ad refers to its value assuming adiabatic activation (i.e., 
e = 0). Entrainment effects on X are also expressed in rel- 
ative terms as 

AX rel = ^% (6) 

Xid 

3.1. CDNC, Q aut , and IF Under Adiabatic Activation 

[ 21 ] The global distribution of CDNC, IF, and Q aut cal- 
culated assuming adiabatic activation (i.e., e = 0) are pre- 
sented in Figure 2 and summarized in Table 2. For the 
meteorological fields used in this study, high CDNC is 
found near heavily polluted areas of the Northern Hemi- 
sphere (NH) (i.e., China, western Europe, eastern coast of 
the U.S. and the Gulf of Mexico), consistent with the high 
level of aerosol associated with industrialized regions, and, 
is the lowest over the polar regions. The influence of long- 
range transport of aerosol is reflected in the higher CDNC in 


the Northern Hemisphere than in the Southern Hemisphere. 
In the latter, high CDNC is associated with biomass burning 
emissions over the continents, and, enhanced sulfate pro- 
duction from DMS oxidation over the oceans. 

[ 22 ] Usage of DAO results in the highest global average 
CDNC (~102 cm 3 ), whereas the lowest is found with 
FVGCM (-95 cm -3 ). Table 2 presents the standard error for 
the comparison, calculated as the mean annual average of 
the standard deviation between the three meteorological 
fields. The large variation in CDNC between the three 
meteorologies (~21 cm -3 ) indicate that the difference 
between global means may not be statistically significant; 
local differences, however, may be much higher. Among the 
meteorological fields used, DAO has the strongest low-level 
poleward transport toward the NH and the weakest wet 
scavenging, resulting in larger aerosol load in the NH than 
with the other fields [Liu et al., 2007]. Similarly, runs with 
FVGCM exhibit the largest CDNC in the Southern Hemi- 
sphere due to larger DMS emissions and stronger poleward 
transport in FVCGM than in the other fields [Liu et al., 
2007], 

[ 23 ] The variability in CDNC and LWC among meteo- 
rological fields results in ~1 ± 0.37 pm variation in R e ff. The 
spread about the mean is, however, much lower than for 
CDNC (Table 2) and the differences are significant to a 90% 
level (given by a student t test). R ef f is typically 1-2 pm 
larger over the ocean than over the continents, which reflects 
the larger LWC (and smaller CDNC) in marine than in 
continental regions [Seinfeld and Pandis, 1998]. These va- 
lues are, however, about 1-2 pm smaller than typically 
observed in satellite retrievals [Bennartz, 2007; Chang and 
Li, 2002] (e.g., Figure 9, section 3.5) although in better 
agreement with ground-based measurements [Dong et al., 
1997], suggesting that while adiabatic activation is reflec- 
tive of cloud-base CDNC, it may lead to overestimation of 
column-average and cloud-top CDNC (hence under- 
estimating satellite-retrieved R e g) [Chang and Li, 2002], 
Local variation in CDNC and LWC results in over 4 orders 
of magnitude variation in Q aut (Figure 2), being about 
1CT 9 s 1 over the southern oceans and as low as 1CT 12 s 1 
in North Africa and East Asia, and higher in the SH than in 
the NH. There is little variation in global mean Q aut among 
meteorological fields; local differences though can be high, 
particularly in the tropics, where DAO and GISS result in a 
factor of 10 higher Q aut than for FVGCM. 

[ 24 ] The global mean IF ranges between -1.28 W nT 2 
(DAO) and -1.75 W nT 2 (FVGCM). As no feedbacks 
between CDNC and LWC are considered in the GMI model, 
the magnitude of the IF predicted is slightly higher than 
current estimates of about -0.7 ± 0.5 W nT 2 [Quaas et al., 
2009]. Thus, IF in Figure 2 represents the first-order 
response (i.e., without feedbacks of CDNC on LWC) of 
shortwave radiative flux to anthropogenic aerosol emissions. 
The largest values of IF are found over SE Asia, western 
Europe and eastern United States consistent with areas of 
high anthropogenic sulfur emissions. High IF (although 
smaller than in the NH) are also found in industrialized 
centers as well as in biomass burning regions of the SH. All 
meteorological fields result in larger IF over the continents 
than over the oceans; long-range transport of pollution 
plumes, however, leads to appreciable IF far into oceanic 
regions, particularly for the DAO and FVGCM fields. The 
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Figure 2. Cloud droplet number concentration (CDNC), autoconversion rate Q aut , and indirect forcing 
(IF) under adiabatic activation (i.e., e = 0). 


runs presented in Figure 2 represent the control case used to 
study entrainment effects on droplet activation. 

3.2. CDNC Under Diabatic Activation 

[ 25 ] Considering entrainment leads to lower CDNC than in 
simulations assuming adiabatic droplet activation (Figure 3 
and Table 3). The global annual mean ACDNC re / lies 
between -10% and -12% for e = 0.4e c , and, from -17% to 
-21% for e = 0.6e c . ACDNC afc (equation (4)) ranges 
between -13 and -15 cnT 3 for e = 0.4e c , and, from -23 to 
-25 cnT 3 for e = 0.6e c ; for the latter the minimum 
ACDNC a/ <,. s is close to -350 cnT 3 (not shown). 

[ 26 ] Consistent with parcel model simulations, areas 
with high CCN concentrations (i.e., high adiabatic CDNC, 
Figure 2), associated with regions of aerosol sources, tend 


to display large absolute changes in CDNC from 
entrainment effects (where ACDNC a /„ is generally around 
-200 cnT 3 ; Figure 3). The largest ACDNC,.,,/ (which is most 
relevant for changes in SW cloud forcing) is, however, found 
in the tropics, downwind of large emission sources, and in 
South America and North Africa (Figure 3). Thus, clouds 
with moderate CDNC (100-300 cm ) are relatively more 
sensitive to entrainment than clouds fonned in very polluted 
regions. 

[ 27 ] The high sensitivity of CDNC to entrainment in 
moderately polluted regions is explained by the sensitivity 
of CDNC to ,v max , which is equal to the slope of the CCN 
spectrum (i.e., the function relating the number of CCN 
with s) at ,v max . CCN spectra for single modal aerosol 
tend to have a sigmoidal shape [e.g., Twomey, 1977; 
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Figure 3. Annual mean ACDNC„,/ for (top) e = QAe, and (bottom) e = 0.6e c . Simulations were carried 
out using (left) DAO, (middle) GISS, and (right) FVGCM meteorology. 
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Medina et al., 2007], the supersaturation range where 
roughly half of the particles act as CCN exhibits the highest 
sensitivity to s (i.e., dC £ N , hence dc ™ c is high). For multi- 
modal aerosol the region of maximum sensitivity to s is 
usually located at lower activation fraction, as the coarse and 
accumulation modes tend to activate at lower s than the 
nucleation mode (Figure 4). For high supersaturations, a 
significant fraction of particles activate so dc £ N — * 0. At low 
supersaturations, the activation fraction and dC £ N are small. 
Thus, if high ,S’ max is reached CDNC is not very sensitive to 
small variations in s max . Figure 4 shows that for the con- 
tinental and marine aerosol distributions of Table 1 the 
region of maximum dc £ N corresponds to an activated 
fraction around 0.1. A second mode of dc £ N is also evident 
around s = 2% about 0.6) which is more prominent 
for the marine aerosol. 

[ 28 ] Since high s max (hence high activation fraction) is 
characteristic of very clean marine environments, 
ACDNC re ; and ACDNC a6x are low there (Figure 5). For 
higher aerosol loads (i.e., near the tropics), CDNC is sen- 
sitive to variations in s max , as the CCN spectrum is steep 


(fraction of particles acting as CCN between 0.2 and 0.5). 
However, for polluted regions (e.g., regions of high sulfur 
emission in the NH), ,s max is low and the CCN spectrum is 
somewhat insensitive to changes in j max (fraction of particles 
acting as CCN below 0.1, Figure 5). ACDNC re / remains 
relatively high at lower activation fraction than ACDNC„/ )S 
as low adiabatic CDNC tends to increase the sensitivity of 
ACDNC re/ to CDNC (equation (4)). 

[ 29 ] Among the meteorological fields used, ACDNC re ; is 
the largest for FVGCM and the lowest for GISS (Figure 3 
and Table 3), reflecting differences in aerosol concentra- 
tion between meteorological fields [Liu et al., 2007]. 
Simulations using DAO and FVGCM result in moderate 
CDNC concentrations (100-300 cm -3 ), hence significant 
ACDNC re / over large regions of the NH. A similar 
spatial variation is found using GISS (Figure 3), although 
ACDNC rf ,/ is relatively lower. Compared to simulations 
with DAO and GISS, the usage of FVGCM fields results 
in higher values of ACDNC,. e / over the southern Atlantic 
Ocean. Figure 5 shows significant scatter of ACDNC a * s 
and ACDNC„,/ for all meteorological fields, reflecting 


Table 3. Global Mean (Maximum) Deviation From Results of Adiabatic Simulation 


Field 


e = 0.4e o 



e = 0.6e c 


e = 2 x 10 3 m 1 
GISS 

e = 3 x 10 3 m 1 
GISS 

DAO 

GISS 

FVGCM 

DAO 

GISS 

FVGCM 

ACDNC„* S (cm -3 ) 

— 14.7(— 151.1) 

—13.4(— 196.8) 

— 1 4.7(— 1 89.8) 

— 25.8( — 26 1 .2) 

-23.5(-345.0) 

-25.8(-333.2) 

-19.5(-210.2) 

— 32.6(— 363. 1 ) 

ACDNC„; (%) 

— 10.4(— 25.4) 

— 1 0.0( — 26.0) 

— 12.1(— 25.8) 

— 1 8.4( — 42.5) 

— 17.7(— 43.2) 

—2 1 ,6(— 42.8) 

— 1 3 ,5(— 55.8) 

—2 1 ,0(— 87.0) 

A R-eff.ab* Om) 

0.25(0.59) 

0.22(0.58) 

0.21(0.58) 

0.48(1.19) 

0.43(1.15) 

0.41(1.16) 

0.33(2.05) 

0.62(3.72) 

QauJ Qaut.ad 

1.12(1.63) 

1.09(1.54) 

1.11(1.58) 

1.26(2.53) 

1.20(2.30) 

1.25(2.40) 

1.18(3.94) 

1.70(32.5) 

A IF abs (W m 2 ) 

0.19(1.18) 

0.16(1.22) 

0.17(1.17) 

0.39(2.21) 

0.32(2.47) 

0.35(2.34) 

0.23(2.13) 

0.40(3.45) 

A IF reI (%) 

— 1 9.6( — 90. 1 ) 

-21.6(— 80.1) 

-15.0(-92.5) 

-31.4(— 99.2) 

-40.4(-98.8) 

-29. 4( -98.4) 

— 25. 1 (—80.5) 

-43.1(-99.5) 
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Figure 4. CCN activation spectra for the continental (solid line) and marine (dashed line) aerosol 
distributions of Table 1. 
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Figure 5. (left) ACDNCq^ and (right) ACDNC,,,/ against activated fraction for e = 0.6e c . Results are 
presented for (top) ocean and (bottom) land. 
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Figure 6. Similar to Figure 3, but for 9- awl . 
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the large variability of aerosol concentration in marine 
and continental environments, and the influence of local 
RH and T. 

3.3. Effective Radius and Autoconversion Rate 

[ 30 ] Compared to adiabatic activation, entrainment reduces 
CDNC and yields a larger effective radius. Thus, AR e jj ahs 
reflects the bias that is introduced when adiabatic activation is 
assumed to represent cloud-column CDNC. Table 3 shows 
the mean and maximum global R e ff for the meteorological 
fields and range of e considered in this study. Large areas 
over the tropics (associated with high LWC and moderate 
CDNC) are characterized by AR e ff abs - 1 //in (e.g., Figure 9). 
For all meteorological fields used in this study global mean 
AR e ff :a h s ranges between 0.21 and 0.48 /./m for e = 0.4e c 
and e = 0.6e c , respectively. 

[ 31 ] According to equation (3), lower CDNC implies 
higher Q aut if q t and C/ are constant. This is shown in 
Figure 6, which presents the ratio of Q aut to its “adiabatic” 
value Qa U t, a d (i.e., for e = 0). Neglecting entrainment effects 
on droplet activation would decrease global Q aut around 
10% for e = 0.4e c and around 23% for e = 0.6e c (Table 3); 
locally, Q au1 can be affected by up to a factor of two. Thus, 
including entrainment effects on CDNC increases the 
autoconversion rate, thus reducing the amount of tuning 
required for a realistic climate simulation [Chen et al., 
2010]. However, if the increased entrainment is accompa- 
nied with an increased rate of droplet evaporation, the net 
increase in Q aut may be negligible [e.g., Jiang et al., 2006; 
Wood, 2007]. Although the global mean Q au JQaut.ad is sim- 
ilar for all meteorological fields used, the regional Q a JQ au t,ad 
is the largest in Central America and northern Africa for 
the DAO simulation whereas for GISS and FVGCM it is 
the largest in the tropical Atlantic and Pacific oceans. 


This is explained by the slightly larger ACDNC,.,,/ in 
marine environments for simulations with the FVGCM 
and GISS fields, compared to the DAO simulations 
(Figure 3). 

3.4. Indirect Forcing 

[ 32 ] The variation of CNDC and R e g with e suggests that 
cloud radiative properties would respond nonlinearly to 
entrainment, therefore affecting IF [Kim etal., 2008], The low 
slope of the CCN spectrum at high .? max (Figure 4) implies that 
CDNC is less sensitive to entrainment during preindustrial 
conditions (because low aerosol loads mean a higher activa- 
tion ratio than for present day). Simulations support this; 
ACDNC„,/ in the preindustrial simulation is lower than for 
present-day simulations, reaching global means around -4% 
and -8% for e = 0.4e c and e = 0.6e c , respectively. Thus, the 
magnitude of IF computed for diabatic activation would be 
lower than when adiabaticity is assumed. 

[ 33 ] Figure 7 presents the annual average A IF abs for all 
meteorological fields of this study. Global mean A IF ahs 
ranges from -0.17 W m 2 to -0.39 W m 2 (Table 3). A//').,./ 
may reach values around -80% in areas of the tropics with 
very low adiabatic IF. In terms of the absolute reduction in 
IF, AIF abs is the strongest (-2.4 W nT 2 ) in regions with 
moderate adiabatic IF (-5 W m 2 ). In pristine regions, 
A IF reI tends to be very sensitive to changes in IF as adia- 
batic IF is very low (equation (4) and Figure 2). These zones 
coincide with low adiabatic CDNC, and AIF rel is the largest 
at very low CDNC, consistent with the high susceptibility of 
cloud albedo to CDNC in clean clouds [Twomey, 1991]. 
However, AIF re/ decreases below -10% for CDNC greater 
than 300 cm 3 (not shown). The largest A IF abs is found at 
moderate CCN concentrations where ACDNC re / is largest 
( — 40%). Because of this, marine environments with mod- 
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Figure 7. Similar to Figure 3, but for AIF ahs . 


erate CCN concentration display a larger A lF ahs (maximum 
around 2.4 W in 2 ) compared to over land where the max- 
imum is ~1.5 W m 2 ( — 35%). A IF abs and A JF reI are small 
at large CDNC because the albedo of very polluted clouds is 
insensitive to variations in CDNC \Twomey, 1991]. 

[ 34 ] The largest global A IF abs (Figure 7) is found for the 
simulations using the DAO field, consistent with the highest 
A R e fiabs of all three fields (Table 3). Global A IF ahs for 
GISS is about 0.03-0.05 W m 2 below DAO values. 
Although the largest AlF abs is located in the tropics, using 
DAO fields tends to concentrate A IF abs near the NH con- 
tinents, whereas for simulations with GISS and FVGCM 
fields, A IF ahs is largest in the tropical oceans. 

3.5. Sensitivity to Variation in e/e c 

[ 35 ] Diagnosing e implies variability in e/e c (section 2.2), 
as e c is determined by the local thermodynamic conditions 
(i.e., p,T,RH). Figure 8 (top) presents the global distribu- 
tion of e!e c for e = 2 x 1(T 3 nT 1 (Figure 8, left) and e = 3 x 
1 0 m 1 (Figure 8, right), using the GISS meteorological 
field; global average e/e c is 0.34 and 0.48, respectively. In 
the tropics, e/e c is considerably above the global average, 
approaching 0.9 for e = 3 x ] 0 3 m ', explained by the 
high T of these regions, resulting in lower e c than in the 
midlatitudes (equation (2)). Thus CDNC in the tropics may 
be more susceptible to neglecting entrainment effects on 
droplet activation than implied in the simulations with 
constant ele c . This is reflected in ACDNC, e / (Figure 8) 
which ranges between -13% and -20%, slightly above 
those obtained using fixed ele c (Figure 3 and Table 3). 
However, ACDNC re / is much larger around the tropics, 
being as much as -70%. This is translated into a mean 
A R e ff,abs about 0.1 pm larger than in the simulations 


assuming constant e/e c (Table 3). In the tropics, A R e ff abs 
can be as high as 3.0 pm (Figure 8). Similarly, global 
mean Q a JQaut,ad is between 1.2 and 1.7, but can be higher 
(particularly in Central America and southern Africa). The 
pattern in IF is also substantially affected. A lF abs (0.23 to 
0.40 W m 2 , Table 3) is globally about 0.1 W m 2 greater 
than assuming a prescribed ele c . This increase results 
mainly from A IF abs in the tropics, which can be as high as 
3.5 W m 2 (Figure 8). 

[36] The reduction in CDNC and increase in R e ff leeward 
of the continents is particularly important near the west 
coastal regions of North America, North Africa, and India 
(also evident in Figure 3) consistent with the CNDC and R e ff 
patterns typically observed in satellite retrievals [e.g., 
Bennartz 2007]. This is further analyzed in Figure 9 where 
the annual mean R e ff using the GISS meteorological field is 
compared against the average R e ff from MODIS retrievals 
for the years 2000-2006 [King et ah, 2006], Compared 
to satellite retrievals, GMI underestimates R^j by about 
3-4 pm which is explained by the higher sensitivity of 
the retrieval to cloud top than to cloud base properties 
[Chang and Li, 2002], Using diabatic activation reduces 
the difference between GMI and MODIS by about 1 pm. 
Qualitatively, diabatic activation results in a global pat- 
tern of R e jf in better agreement with MODIS, than when 
adiabaticity is assumed. This is particularly noticeable in 
the remote ocean regions of the tropics (0 to 30°S) for e = 
3 x 1 0 3 m ', where regions of large R,.ff across the 
tropical oceans are much better captured by GMI. The 
contrast in R e ff between land and ocean in the midlatitudes 
is also better represented using diabatic activation (e = 
2 x 1CT 3 ), particularly off the West Coast of North 
America, and in the east coast of Africa and South America. 
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Figure 8. Annual global distribution of e/e c , ACDNC re i, AR e ff,abs, q ““‘ an( J AIF ahs for e 
(left) m _I and (right) e = 3 x 10 -3 m _1 , using GISS meteorology. 
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The characteristic e for diabatic droplet activation may 
depend on the local wind dynamics [Siebesma and Cuijpers, 
1995]. Figure 9, however, suggests that its variation may be 
small over large areas of the globe. 

4. Discussion and Conclusions 

[ 37 ] The global distribution of CDNC, IF, R e /j, and Q aut 
for diabatic droplet activation was studied with the NASA 
GMI model using three different meteorological fields. In 


general, considering diabatic activation resulted in lower 
CDNC than when adiabatic activation was assumed. The 
largest variability in CDNC (up to -75%) was found in 
the tropics and in zones of moderate CCN concentration. 
The redistribution of LWC into a lower diabatic CDNC 
resulted in global mean ./^between 0.2-0. 5 /tm (and locally 
up to 3.5 //m over the tropical oceans) higher than for adi- 
abatic activation. Similarly, the global mean Q aut is increased 
by a factor of 1.1 to 1.7, and as high as a factor of 4 in the 
tropics. It was found that by assuming adiabatic activation, 
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Figure 9. Annual global distribution of R^j ijnw) calculated using the GISS meteorology. Also shown 
are A 1 ,,// retrievals from the NASA MODIS platform averaged over the years 2000 to 2006 [ Meskhidze 
et al., 2007; King et al., 2006]. 


IF may be overestimated between 15 and 40%. Entrainment 
effects on droplet activation and indirect forcing were more 
significant using the DAO meteorological field as AR e ff abs 
was typically higher for these simulations, which is linked to 
the higher LWC in tropical clouds predicted by these field 
than with FVGCM and GISS. 

[38] It was shown that areas with moderate CCN con- 
centration (i.e., activated fraction between 0.2 and 0.5) are 
most prone to entrainment effects during activation (and 
therefore may produce a larger bias in R e ff when adiabatic 
CDNC is assumed) as the CCN spectrum is at its steepest, 
and CDNC most sensitive to variations in s m ax . Using 
adiabaticity as a proxy for e/e c assures that our approach is 
consistent with current observations of the entrainment 
effect on LWC [e.g., Lu et al., 2008; Kim et al., 2008; 
Morales et al., submitted manuscript, 2010]. However, since 
e/e c may vary throughout the globe, we considered simu- 
lations for which e was constrained using published LES 
simulations and field observations. Not surprisingly, IF and 
CDNC can be sensitive to the value of e used. Increasing e 
from e = 2 x 10 3 m 1 to e = 3 x 10 3 m 1 would decrease 
global mean ACDNC,.^ from -13% and -20% and increase 
the global mean AIF ahs from 0.23 to 0.40 W m 2 . Thus, 
neglecting entrainment in droplet activation parameteriza- 
tions would add up to about 36% uncertainty to estimations of 
the aerosol indirect effect. As e can vary between e = 0.5 x 


1(T 3 m _1 and e = 4 x 10 m _1 [Siebesma and Cuijpers, 

1995; McCarthy, 1974; Raga et al., 1990; Gregory, 2001; 
Neggers et al., 2003], variations in the parameter could impact 
IF assessments (e.g., e = 4 x 10~ 3 results in about 50% lower 
IF than for adiabatic activation). Hence entrainment effects 
on droplet activation (particularly in the tropics) may need 
to be considered in GCM studies. 

[ 39 ] The contrast in CDNC and R e ff between continental 
and marine environments found in satellite retrievals [e.g., 
Bennartz 2007] is much better captured when diabatic 
(Figures 3 and 9) rather than adiabatic (Figure 2) activation 
is assumed. Using diabatic activation substantially improved 
the prediction of R e ff when compared to satellite retrievals, 
particularly in marine environments where the effect of 
entrainment on CDNC was the strongest (Figure 9). It is 
remarkable that the simple entrainment correction presented 
in equation (1) resulted in such an important improvement 
on the global pattern of R e ff predicted by GMI. 

[ 40 ] Although this study is focused on the study of the 
first indirect effect, the impact of diabatic activation on Q aut 
suggests that the cloud lifetime can also be affected. The 
larger R e ff for diabatic than for adiabatic activation would in 
principle lead to larger precipitation rate and shorter cloud 
lifetime. However, this would feedback on LWC limiting 
Q aU f Partial droplet evaporation from entrainment (not 
considered in this study) is also possible and would 
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affect the droplet size distribution and Q aut [Brenguier 
and Chaumat, 2001]. Including dynamical and hydrologi- 
cal feedbacks from changes in the radiative balance and 
autoconversion can dampen (or magnify) the magnitude of 
CDNC responses to entrainment [Stevens and Feingold, 
2009], and is required to assess the effect of diabatic acti- 
vation on cloud lifetime. Still, the results of this study 
suggests that diabatic effects on droplet formation may be an 
important omission in global model studies of the aerosol 
indirect effect. 

[ 41 ] In this study aerosol distributions are prescribed and 
scaled to aerosol constituent mass. Although this approach 
is known to introduce uncertainties in the predicted CCN 
spectrum [Meskhidze et al., 2007], particularly far from 
sources where aging and mixing can be significant, the 
conclusions of the present study may not be as sensitive to 
them, as we are studying the relative impact of assuming 
adiabatic/diabatic activation on CDNC, R e jj, and IF. We 
have also favored the possibility of running a large number 
of cases over using a high-resolution GCM. Using a finer 
resolution would not change the conclusions of this study, 
however, may modify the absolute values of CDNC, R e jj, 
Q aut and IF. All these issues are important and require the 
application of a fully coupled climate model with explicit 
aerosol microphysics, and will be addressed in future 
studies. Still, the magnitude and sign of ACDNC,,/„ and 
A Reffabs found in this study point in the right direction as 
R e ff is typically underestimated and CDNC overestimated 
by GCMs that use adiabatic CDNC values for forcing 
calculations. Thus, considering diabatic activation can 
provide a first-order correction to CDNC calculations that 
lead to more realistic distributions in predicted cloud 
microphysical properties and in the estimation of the 
aerosol indirect effect. 

Notation 

Qaut, Qaut, ad autoconversion rate computed using diabatic 
and adiabatic activation, respectively, 
abs, rel subscripts indicating the absolute and relative 
change, respectively, from adiabatic (e = 0) to 
diabatic activation. 

CCN cloud condensation nuclei. 

C/ liquid water cloud fraction (stratiform). 

CN condensation nuclei. 
c p specific heat capacity of air. 

e, e e ff effective fractional entrainment rate. 
e c critical entrainment rate, 
g acceleration of gravity. 

IF indirect forcing. 

M w ,M a molar masses of water and air, respectively. 
p Pressure. 

p\T) saturation vapor pressure of water. 
qi liquid water mixing ratio (stratiform). 

R universal gas constant. 

Reff cloud droplet effective radius. 

RH ambient relative humidity. 
s water vapor supersaturation. 

Smax maximum water vapor supersaturation in a 
cloud parcel. 

T cloud parcel temperature. 
t time. 


T' ambient temperature. 

V updraft velocity. 

W liquid water mixing ratio. 

gAH v M„. _ gMa_ 

^ c p RT- RT ' 

pM a , M W AH 1 
'7 p , (T)M w CpRT 1 • 

A H v enthalpy of evaporation for water. 
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